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Abstract. We systematically study the effects of collisions on the overall dynamical evolution of 
dense star clusters using Monte Carlo simulations over many relaxation times. We derive many 
, observable properties of these clusters, including their core radii and the radial distribution of 
collision products. We also study different aspects of collisions in a cluster taking into account the 
shorter lifetimes of more massive stars, which has not been studied in detail before. Depending 
on the lifetimes of the significantly more massive collision products, observable properties of 
the cluster can be modified qualitatively; for example, even without binaries, core collapse can 
sometimes be avoided simply because of stellar collisions. 
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P 1. Introduction 

c/3 In dense stellar systems like massive young clusters, galactic centers, and old globular 

clusters (GC), the high densities of stars can give rise to many direct, single-single 
physical collisions, in addition to collisions mediated by dynamical interactions of binaries 
( [Fregeau et al . 2004) . In these systems stellar evolution can also be significantly modified 
through physical collisions. For example, in young star clusters, collisions can give rise 
to runaway growth of merger products, producing many exotic stellar populations, such 
as intermediate- mass black holes (Giirk an et al. 20^61 ITrenti 2006p . The cores of typical 
old Galactic GCs can also attain high enough densities so that most core stars undergo 
collisions during their lifetime ( [Hills fc Day 1976} [Fregeau et al. 2004p . These collisions 
^-H not only change the evolution of individual stars, but the increased stellar masses also 
C — , change the overall GC properties like the core radius (re) and the half-mass radius {rh)- 
■ Moreover, at least some of the observed exotic stellar populations in dense clusters, 
j>! ' like blue straggler stars (BSS) (Feka du et al. 2007HSiirs et al. 200lllSills et al. 1997p . and 
' compact binaries like ultracompact X-ray binaries (UCXB), are likely created through 
[ collisions ([Lombardi et al. 2006p . 

' A dense cluster of stars naturally evolves towards eventual core collapse through re- 
laxation. At their present ages, most Galactic GCs are expected to have collapsed cores. 
However, observations show that the measured values of Vc/rh for most Galactic GCs 
are higher than predicted by theoretical models (Vesperini fc Chernoff 1996| . Many sce- 
narios have been proposed to explain this apparent discrepancy. For example, the core 
can be supported against deep collapse by dynamically extracting the binding energy of 
hard primordial binaries ([Trenti et al. 20071 [Fregeau fc Rasio 2007[ ). However, it is hard 
to explain most of the high observed values of rc/vh in Galactic GCs purely through 
this "binary burning" process. Other mechanisms to halt core collapse like ejection of 
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Table 1. Initial Conditions of Simulated GCs 





Model 


IMF 


N 


f^binary 


Tc (pc) 


Virial radius (pc) 




Cases 1,2 


Plummer 


Single-Mass 


10" 





0.3 


0.85 


10" 


Case 3 


King, Wo = a 


Kroupa (0.1 - 2.0Mq) 


10" 





0.87 


2.89 


1.7 X 10" 



stellar mass BHs from the core (Mackey et al. 2007) or stellar captures by a central 
intermediate-mass black hole (jTrenti et al. 20071 ITrenti 2006|) have also been discussed 
at this meeting (see contributions by Mackey and Trenti in this volume) . It has also been 
suggested that rdrh could simply keep increasing over long timescales in clusters having 
fairly long relaxation times via ongoing mass segregation (jMerritt et al. 2004j) . 

Here we study stellar collisions as a possible mechanism for supporting clusters against 
core collapse. In a regime where collisions are important, they can produce many stars 
significantly more massive than those in the background population. Thus the subsequent 
evolution of collision products will be much faster than for normal stars in the cluster. 
Although the stellar evolution and observable properties of collision products have been 
extensively studied before (S ills et al. 20011 fSills et al. 1997|) . the feedback effects of col- 
lisions on the overall dynamical evolution of GCs has received less attention. The shorter 
evolution timescales of massive collision products may support the core against collapse 
even without any primordial binaries or other mechanisms for energy production, simply 
via indirect heating through stellar evolution mass loss ( [Goodman fc Hernquist 1991[ 
IGoodman fc Hut 19891 iLee 19871) . 

Using iV-body simulations, we have begun studying numerically how collisions and the 
subsequent evolution of collision products can alter the overall properties of GCs. We use 
the Northwestern group's Henon-type Monte-carlo code, which provides a detailed, star- 
by-star representation of clusters with up to TV '--^ 10^ — 10'' stars ([Fregeau &: Rasio 2007 



IJoshi et al. 200l|IJoshi et al. 2006|) . In §2 we will first present two simple limiting cases, 
bracketing reality and illustrating the dramatic changes in global cluster properties de- 
pending on how the evolution of collision products is treated in the models. We also 
present the evolution of a more realistic GC model with a conventional "rejuvenation" 
prescription for determining the lifetimes of collision products. We discuss the implica- 
tions of our study and planned future work in §3. 



2. Results 

All simulations shown here use N — 10^ stars initially. We assume totally conservative 
collisions with the interaction cross section given by the usual "sticky sphere" approx- 
imation. There is no hydrodynamic mass loss during the collisions. However, mass loss 
occurs through normal evolution of the stars via (instantaneous) compact object for- 
mation at the end of the main sequence (MS) lifetime. The collision products can be 
significantly more massive and hence evolve faster, enhancing the overall stellar mass 
loss rate. Determining theoretically the true MS lifetime or the effective age of a collision 
product is a hard task. These quantities depend sensitively on the amount of hydrogen 
mixed into the core of the collision product as a result of the collision. This in turn can 
be highly variable because the mixing depends on the details of the collision kinematics. 
To avoid these uncertainties, in our study we use three very simple prescriptions for de- 
termining the effective ages and the remaining MS lifetimes of collision products. Case 
1: Infinite lifetime for all stars, including the collision products. Case 2: Zero lifetime 
for the collision products. Other stars have infinite lifetime, unless they collide. Case 3: 
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Figure 1. a) Evolution of r^/rh- Solid line: no collisions. Dotted line: case 1, where all stars 
have infinite lifetime (including collision products). Dashed line: case 2, where collision products 
have zero lifetime. In case 2 collapse is avoided through mass loss due to evolution of the collision 
products, b) Evolution of E (top panel) and the core radius (bottom panel). The steepest slope 
in the E evolution curve corresponds to core collapse. Onset of collapse can also be seen directly 
in the rc evolution plot, c) Evolution of E (top panel), evolution of rc in physical units (middle 
panel), and Evolution of the average mass contained in three Lagrange radii (0.01,0.1,0.9) of 
the GC (bottom panel) for case 3. d) Evolution of rc/r^ for case 3. e) Histogram of rc/r^ of all 
Galactic GCs with known Vc and values. 

A more realistic "rejuvenation" prescription for determining the age and MS lifetime of 
the collision products: we assume that stars burn H linearly on the MS so the effective 
age of a collision product is uniquely determined by the total amount of H coming from 
the progenitors. The age and the MS lifetime of each collision product is determined by 
the following simple equation, based on Hurley et al. (2002) [ 

tMS fhMi _ t2M2' 



t = 



(2.1) 



M \tMSl tMS2 , 

Here, U, Mi, and Imsi are the age, mass, and the MS lifetime of parent star i, while t, 
M, and tMS are for the collision product. Note that the coefficient on the right hand side 

We use Eq.El] so that 



of Eq. l2.1l is different from the one used in Hurley et al. (2002) 



two zero age main sequence (ZAMS) stars collide to produce a ZAMS star, whereas two 
stars close to their MS turn-off collide to produce another star close to its turn-off. 

2.1. Limiting cases: Single mass Plummer model (cases 1 & 2) 

The initial cluster properties are listed in Table[T] Fig.[T^ shows the evolution of rc/vh 
for the two limiting cases as well as the case without collisions for comparison purposes. 
Clearly, conservative collisions make the cluster more bound manifested by the faster 
collapse of the core than the no collision case in case 1. In the other limiting case, case 
2, on the onset of collapse the densities reach a very high value increasing the collision 
rate. The mass loss after the collisions can stop the collapse and reach a steady core 
radius with a low rate of further collisions. This mechanism can be better illustrated by 
Fig-Hb- We define energy E in such a way that it remains constant unless mass is lost 
from the cluster. The onset of collapse can be seen in the simultaneous decrease of rc 
and the increased slope of the E evolution curve. 
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2.2. More realistic case: rejuvenation- based prescription; King sphere, wq — 6 (Case 3) 

The initial cluster properties are listed in Table[T] Here we first evolve the cluster with 
only stellar evolution for 10^ years, so that the very massive stars evolve and disappear 
from the system. We then switch on dynamics. 

Till a little more than 1 Gyr E remains constant, indicating no collisions (Fig.[TJ;). 
Sometime between 1 & 2 Gyr, the core starts to collapse (Fig.[T];) reducing the core 
radius and increasing the average mass in the core region. Thus collision rate increases 
significantly increasing rate of mass loss. This is manifested in the sudden steep increase 
in E. Thus through mass loss core collapse is avoided and rc/r^ approaches a steady 
value (Fig[TJi). Fig[TJl,e shows that, even with the simple rejuvenation prescription, the 
final value of rc/r^ compares well with the observed values for most Galactic GCs. The 
slow continued increase of E at later times comes as normal stars as well as some low 
mass collision products go off the MS and disappear. 

Fig. [2^ shows the positions of the collision products that are still in their MS life at 
14 Gyr. Most of these collision products are contained within the core. BSS candidates 
are most easily identified above ~ 2 times the turn-off mass of the cluster. The more 
massive subset of collision products in our simulated cluster (M ^ I.QMq) are mostly 
found inside the core (Fig.[2j:) (apart from only 2 just outside the core). 



3. Discussion 

Although determination of stellar evolution after a collision is a subject of continuing 
research, the effects of collisions on cluster dynamics have not been studied in detail 
previously. This is a first report of an ongoing systematic study on this effect. Using 
typical initial conditions for old GCs and simple assumptions for rejuvenation we show 
that collisions between single stars not only alter the stellar properties and produce exotic 
stellar populations like some of the BSS, they also affect the overall GC properties. Most 
importantly, collisions can support the core of a cluster against collapse in typical old 
Galactic GCs. Even with our simple but reasonable assumptions, the values of rc/rh 
obtained for our simulated GCs compare well with the observed rc/rh values for Galactic 
GCs. Furthermore, we obtain a population of BSS candidates contained within the core, 
also consistent with observations (Leigh et al. 2007). Note that we do not find any BSSs 
well outside the core, consistent with the current understanding that those BSSs arc most 
likely formed via primordial binary mergers ( [Mapelli et al. 20"06l ) . 

We have adopted many extreme simplifications for this first look at the problem. For 
example, at this stage of our study, once a star evolves off the MS, it is removed from 
the simulation, leaving no remnant. Since remnants are normally only a few percent of 
the total progenitor mass, we expect that this approximation will not affect the over- 
all GC properties significantly, so far as the increase in energy is concerned from mass 
loss. However, some remnants from very massive stars at a very young age can remain 
in the cluster and sink into the core through dynamical friction. Dynamical interac- 
tions, including collisions, of these massive remnants can also alter the core properties 
of the GCs in certain regimes ( |Mackey et al. 2007| ITrenti et al. 20071 [^nti 2006p . An- 
other possibly important effect left out of this study for now is the role of primordial 
binaries ( [Fregeau k. Rasio 2007[ [Fregeau et al. 2004] ). On the one hand, the presence 
of primordial binaries will increase the rate of collisions through resonant encounters 
( [Fregeau et al. 2004p , on the other hand, binaries will provide further support of a clus- 
ter against collapse and hence may prevent the core from reaching high enough densities 
for significant collisions. 
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Figure 2. a) Position vs mass scatter plot for all collision products still at their MS life at 
14 Gyr. Two times the present age turn-off mass (1.6Mq) is shown as a horizontal solid line 
to guide the eye. The position of the core is shown with the vertical dotted line, b) Histogram 
showing the positions of the same population as in (a), cj Histogram showing the positions of 
the collision products still in their MS life at present age having masses ^ 1.6Mq. 
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